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INTRODUCTION 


The  research  on  sensitivity  coefficients  on  ballistics  has  been  discussed 
in  two  previous  works.  The  first  paper^  matches  the  sensitivity  coefficient 
of  exterior  ballistics  to  that  of  interior  ballistics,  while  the  second  paper? 
presents  the  sensitivity  coefficient  of  exterior  ballistics  with  velocity 
square  damping.  This  report  continues  the  work  of  the  latter?  and  extends  the 
analytic  study  in  order  to  give  numerical  solutions  of  the  problem. 

The  design  of  a  gun  Involves  numerous  parameters.  These  parameters 
.should  be  in  such  a  combination  that  the  best  first  round  accuracy  is  given. 
While  the  shell  leaves  the  gun  it  has  perturbations  for  the  muzzle  elevation 
angle  and  the  muzzle  velocity.  The  ratio  of  the  two  is  the  sensitivity  coef¬ 
ficient  of  the  interior  and  the  exterior  ballistics.  It  is  desired  to  compen¬ 
sate  for  errors  due  to  uncertain  changes  of  muzzle  velocity  by  the  automatic 
response  of  the  muzzle  elevation  angle  within  the  gun  system.  With  a  correct 
design  this  compensation  can  be  made  by  matching  the  exterior  ballistics  to 
the  interior  ballistics  through  the  analysis  of  gun  dynamics.  This  process  is 
called  passive  control  since  there  is  no  external  measurement  Involved  nor 
Instrumentation  needed  for  control.  This  general  problem  can  be  formulated  by 
first  investigating  the  sensitivity  coefficients  for  exterior  ballistics  with 
velocity  square  damping. 

Istien,  C.  N.,  "On  the  Sensitivity  Coefficient  of  Exterior  Ballistics  and  Its 
Potential  Matching  to  Interior  Ballistics  Sensitivity,”  Proceedings  of  the 
Second  US  Army  Symposium  on  Gun  Dynamics  at  the  Institute  on  Man  and  Science, 
Rensselaervllle,  NY,  September  1978,  sponsored  by  USA  ARRADCOM. 

?Shen,  C.  N.,  "Sensitivity  Coefficient  of  Exterior  Bal.  1 sties  With  Velocity 
Square  Damping,"  Transactions  of  the  Twenty-Fifth  Conference  of  Army 
Mathematicians,  ARO  Report  80-1,  January  1980,  pp.  267-282. 
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PREVIOUS  WORK 


In  the  second  paper  the  principal  equation  of  exterior  ballistics  has  a 
drag  term  which,  in  this  case,  is  proportional  to  the  square  of  the  velocity 
In  the  tangential  direction  of  the  projectile.  The  sensitivity  coefficient  Is 
expressed  as  the  ratio  of  the  Initial  elevation  angle  deviation  to  the  initial 
precentage  velocity  deviation.  The  work  Is  to  find  analytically  the  sensitiv¬ 
ity  coefficient  of  the  exterior  ballistics  with  velocity  square  damping  which' 
comes  from  the  nonlinear  air  resistance  for  a  projectile.  This  principal 
equation  is  integrated  analytically  in  obtaining  the  solution  for  tangential 
velocity  in  terms  of  the  elevation  angle,  together  with  all  the  necessary  ini¬ 
tial  conditions.  The  horizontal  range  and  the  vertical  range  are  also 
expressed  as  integrals  of  certain  function  of  the  elevation  angles.  In  order 
to  obtain  the  sensitivity  coefficient  it  is  necessary  to  find  the  perturba¬ 
tions  of  the  horizontal  and  vertical  ranges.  This  procedure  is  similar  to 
that  of  evaluating  differentiation  under  the  integral  sign.  The  perturbation 
of  the  ranges  is  the  sum  of  the  perturbations  due  to  the  initial  velocity,  the 
initial  elevation  angle  and  the  Impact  elevation  angle.  By  setting  to  zeroes 
the  range  perturbations  we  can  group  the  coefficients  of  the  perturbations 
into  two  separate  equations.  The  ratio  of  the  perturbations  for  initial 
elevation  angle  to  that  for  initial  velocity  is  the  sensitivity  coefficient 
for  exterior  ballistics  that  we  are  seeking. 


RESULT  OF  PAST  RESEARCH 

The  following  results  are  concluded  as: 

1.  The  principal  equation  of  exterior  ballistics  la  derived  with  the 
trajectory  slope  as  the  Independent  variable. 

2.  The  closed  form  solution  for  the  horizontal  component  of  trajectory  ■ 
velocity  is  determined  for  the  case  of  exterior  ballistics  with  velocity 
square  damping. 

3.  The  nondimenslonal  range  Is  obtained  In  terms  of  an  end  slope 
function  and  a  range  drag  function. 

4.  Variations  of  the  nondimenslonal  range  are  expressed  as  variations  of 
launch  velocity. 

5.  Variations  of  the  range  drag  function  are  in  terms  of  the  variations  ' 
of  the  range  drag  Integral. 

6.  The  range  drag  integral  has  parameters  in  the  integrand  as  well  as 
the  upper  and  lower  limits.  The  variations  of  the  integral  are  found. 

7.  The  partial  derivatives  of  the  range  drag  integrand  are  evaluated. 

8.  The  variational  equation  for  the  range  Is  In  terms  of  elements 
Involving  three  Integrals  as  coefficients  of  three  variational  parameters. 

9.  The  variational  parameters  are  that  of  launch  velocity,  the  launch 
elevation  angle,  and  the  Impact  elevation  angle. 

10.  The  average  of  the  end  slopes  Is  equal  to  the  terrain  slope  times  the 
range  drag  function  minus  the  elevation  drag  function. 

11.  Variations  of  the  nondimenslonal  elevation  are  expressed  as 
variations  of  the  end  slopes  and  the  variations  of  the  drag  function. 


12.  The  variational  equations  for  the  elevation  are  determined  similar  to 
that  for  the  range. 

13.  Eliminating  the  variations  of  impact  slope,  6qj ,  from  the  set  of  two 
variational  equations  gives  the  ratio  of  the  coefficients  of  6v0/v0  and  6q0/ 

(q0"<ii>* 

14.  The  sensitivity  S80/(6v0/v0)  may  be  obtained  by  dividing  this  ratio 
6q0/(5v0/v0)  by  the  quantity  (l+q02). 

However,  numerical  calculation  of  this  problem  was  not  carried  in  that  paper.2 

DYNAMICAL  EQUATIONS  FOR  TRAJECTORIES 

After  several  transformations  as  given  in  the  paper2  and  summarized  in 
Appendix  A,  the  dynamical  equations  are  simplified  as 

du  -  -  u3  (l+q2)l/2dq  (1) 

q 

u2 

dx  -  -  —  dq  (2) 

g 

u2 

dy  - - qdq  (3) 

g 

In  the  above  equation  the  independent  variable  q  is  the  projectile  slope  which 
is  related  to  8  as 

q  -  tan  8  (4) 

and  the  dependent  variable  u  is  related  to  v  and  8  as 

u  -  v  cos  8  (5) 

2Shen,  C.  N. ,  "Sensitivity  Coefficient  of  Exterior  Ballistics  With  Velocity 
Square  Damping,"  Transactions  of  the  Twenty -Fifth  Conference  of  Army 
Mathematicians,  ARO  Report  80-1,  January  1980,  pp.  267-282. 
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-•-SiSili 


where  g  *  the  acceleration  due  to  gravity 


c  -  the  drag  coefficient  of  the  projectile 
v  •  the  velocity  of  the  projectile 
0  ■  the  path  inclination  (elevation  angle) 
x  •  the  horizontal  distance  of  the  projectile 
y  -  the  altitude  or  vertical  distance  of  the  projectile 
The  solution  for  horizontal  component  of  velocity  u  can  be  obtained  by 
Integrating  Eq.  (1)  to  give 


where 


and 


— 2  ■  — o  U"uo2  “  tp(q)  -  P0(q0)]} 

UZ  Uq^  g 


p(q)  -  qd+q2)1/2  +  Mq  +  (1+q2)1/2] 
P0(q0)  ••  qod+qo2)1^2  +  *ntq0  +  d+qo2)1^2! 
UQ2  -  vo2sec"20o  -  v02(l+q02)"1 


q0  -  tan  90 


Finally,  Eq.  (6)  takes  the  form 

„  2 


u2  - 


(1  + 


H(q,q0*v02 »c/g) 


i+qQ2  1  "  H(q»qo»vo2*c/s 


where 


c  vo 


w(q»q0»vo2»c/8)  ■  -  r — ?  "  Po(q0)l 

s  i+q0 


(6) 

(7) 

(8) 

(9) 

(9a) 


(10) 

(11) 


SOLUTION  FOR  NONDIMENSIONAL  RANGE 

In  determining  the  range  x  for  the  trajectory  the  closed  form  solution  of 
u2  in  Eq.  (10)  can  be  substituted  into  Eq.  (2)  to  obtain  the  solution  in 
integral  form  as 


j 


•I  L  . •■"^fP<|WPS|| 


,7r^>  I<,,',,,)  + 

where  x^  ■  range  at  Impact  point 
xQ  ■  range  at  Initial  point 
and  qj  ”  projectile  slope  at  impact  point. 

To  non-dimensionalize  the  range,  Eq.  (12)  ia  divided  by  the  factor  v02/g 
as 

X(xltx0,v0)/A(qi,q0)  -  GxUi ,q0.v02 »c/g)  03) 

where  the  nondimensional  range  is 

X(xi,x0,v0)  -  (xi-x0)g/vo2  (14) 

the  slope  function  is 

q0-qi 

Mq0»qi>  *  05) 

i+q0 

and  the  range  drag  function  due  to  air  resistance  is 

1  Qi  H(q,q0,v02,c/g) 

Gx(qi»qo»vo2»c/8)  -  i - I  - 7 — - - d<l  06) 

q0-qi  q©  1  -  h 


qi  H(q,q0,v02,c/g) 

/  - dq 

q0  1  -  H 


(12) 


VARIATION  OF  THE  NONDIMENSIONAL  RANGE  AND  THE  SLOPE  FUNCTION 

In  order  to  obtain  a  first  round  hit  of  the  target  one  of  the  conditions 
is  that  the  variation  of  the  range  should  be  zero,  i.e.,  from  Eq.  (14) 

6(xi-x0)  -  0  (17) 

We  take  the  perturbation  for  the  nondimensional  range  from  Eq.  (14)  as 

6X  5<*i-*o>  2«v0 

X  xt-x0  v0 


I 


(18) 


fix  2fiv0 

...  .  .... 

The  variation  of  the  slope  function  A  in  Gq.  (15)  becomes 

fiA  5qi  $<lo  2q0fiq0 

A  <lo-<U  lo-ll  l+<lo2 


(19) 


Next,  taking  the  variation  of  Eq.  (13)  and  using  the  expressions  given  in  Gq. 
(18)  and  (19)  we  have 


fiX  fiA 

fiGx 

— 

(20) 

X  A 

fiGg 

26v0  6qt 

$<lo  2q0fiq0 

(21) 

— - — —  +•  — — — — -  - 

”  — — — — —  +  -••••• 

Gx 

vo  90-91 

qo“9i  i+q0 

This  gives  the  variation  of  the  range  fiGg  in  terms  of  the  variations  of  the 
initial  velocity  fivQ>  and  the  variations  of  the  Initial  slope  fiqQ  and  that  of 
the  impact  slope  fiqj. 


THE  SOLUTION  FOR  ELEVATION 

The  differential  equations  for  elevation  was  given  in  Eq.  (3)  and  the 
solution  for  u  is  in  Eq.  (10). 

Substituting  Eq.  (10)  into  Gq.  (3)  gives 


v02  H(q,q0,v02,c/g) 

dy  - - — - x-  11  - - - Iqdq 

S ( l+q02 )  1  *  « 


Integrating  the  above  one  obtains 

v  2 
vo 

yi-y0  -  - 


9i2-q02  Ri  qH 

- - -  [ - +  /  - dq] 

g(l+q02)  2  q0  1-H 


(22) 


(23) 


7 


Rearranging  yields  the  relationship  between  the  range  Y,  the  end  slope 
function  A,  and  the  elevation  drag  function  Gy. 

1  , 

Y<yi,y0tV0)/A(qltq0)  -  -  (q0+qj)  +  Gy^  ,q0,v0%c/g)  (24) 

where  the  nondimensional  elevation  is 

g(yi-Yo> 

Y(yi.y0.v0)  -  — 5 —  (25) 

A  is  given  in  Gq.  (15),  and  the  elevation  drag  function  is 

1  <li  qH(q,q0,v02,c/g)dq 

Gy(qi,q0,v0z,c/g)  - - /  - - -  (26) 

q0-qi  qo  I-® 

TERRAIN  SLOPE  FROM  LAUNCH  POINT  TO  TARGET  POINT 

If  Eq.  (25)  Is  divided  by  Eq.  (14)  with  the  ?ld  of  Eqs.  (16)  and  (26), 
one  obtains 

Yl-Yo  (l/2)(q0+qi)  +  Gy 

- Am- - -  (27) 

*i"*o  -  Gx 

where  m  Is  the  terrain  slope  from  launch  point  to  target  point,  a  constant 
parameter.  Therefore, 

(l/2)(q0+q£)  +  Gy  —  aGj  (28) 

We  use  Eq.  (27)  to  find  the  variational  equation  for  the  elevation. 

Taking  the  variation  of  Eq.  (28)  for  any  given  m,  we  have 

(i/2)(6q0+6q1)  +  «Gy  -  m  6G*  -  0  (29) 

It  Is  noted  that  6Gy  and  dGx  are  related  In  Eq.  (29),  with  also  the  variations 
6q0  and  6q£. 
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VARIATION  OF  THE  RANGE  AND  ELEVATION  DRAG  FUNCTIONS 


Work  was  performed  on  evaluation  of  the  drag  function  In  the  previous 
paper^  which  Is  summarized  In  Appendix  B.  This  gives  the  respective  varia¬ 
tions  of  the  range  and  elevation  drag  functions  as 


where 


26v0 

6qt 

«q0 

vo 

+  ai 

T  “0 

qo*qi  9o-9i 

2Sv0 

■v  hi 

6qt 

' 

Sq0 

u 

vo 

-r  U0 

qo-qi  q0-qi 

1 

av  -  [ - I12] 

do-<li 


(30) 


(31) 


(32) 


Hi  1 

ai  “  t-  - - Ill! 

i-Hi  q0-qi 


(33) 


2q0  1  c  vo  dPo 

a0  -  [  r  1 1 2  + - Ill  +  (")  ~~~  l02l 

i+q02  qo+Hi  8  1+Ho2  dqQ 


(34) 


bv  -  [ - Jl2l 

q0-qi 


(35) 


qiHi  1 

bi  ■  l-  7-- - dlll 

i-Hi  qo'qi 


(36) 


and 


2q0  1  c  vo  dPo 

5o  -  l  77--J  J12  + - J11  +  (-)  ~7~"  T~  J02l 

i+qQ  q0-qi  8  i+q0  d<io 


(37) 


2shen,  C.  N. ,  "Sensitivity  Coefficient  of  Exterior  Ballistics  With  Velocity 
Square  Damping,"  Transactions  of  the  Twenty-Fifth  Conference  of  Army 
Mathematicians,  ARO  Report  80-1,  January  1980,  pp.  267-282. 
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In  turn,  the  integral  In,  In*  and  lQ2*  ottar  terms  are  given  as 


follows. 


and 


qi  H(q,qo,v02,c/g) 

in(qo.qi)  -  /  - t~ - d<» 

q0 

Il2«lo.1l)  •  /’*  “2  “I 

^“2  *■ 

%  •  H(q-qi) 


dpo 

dqo 


20+qo2)1/2 


<31 

(3! 

(4< 

(4 

(4 


Moreover,  the  Integral  Ji2»  ai*H  J02»  an<*  °tHer  terms  are  given  as 

follows. 

qi 

Jll  -  /  qH(q,q0,v02 ,c/g)/ (1-H)  dq 

qo 


(4 


SOLUTION  FOR  SENSITIVITY  COEFFICIENT 


It  1 8  noted  that  both  Eq.  (21)  and  Eq.  (30)  give  the  variation  fiG*  in 
terms  of  variations  6vQ,  6qt ,  and  6q0  and  the  equations  are  Independent  of 
each  other.  Rewriting  these  equations  as: 


26v0 

dqj 

«9o  2q06q0 

SG*  “  Gx  t . 

+ - 

— — — —  +  — — ?] 

(47) 

vo 

9o-9i 

q0-<U  l+9o 

26v0 

6qi 

$q0 

6Gx  -  av - 

+  ai - 

-  +  ao - 

(48) 

vo  qQ-qi  ‘lo-^i 


Subtracting  Eq.  (48)  from  Eq.  (47)  and  then  dividing  by  Gx  gives 

av  26 v0  ai  6q*  aQ  6qQ  2q06q0 

0  -  [-1  -  --] - +  [1  -  --] - +  [-1  -  --1 - +  ~~l 

Gx  vo  Gx  qo"*qi  Gx  qo~qi  1+90 


(49) 


It  is  also  noted  that  Eq.  (31)  expresses  the  variation  6Gy  in  terms  of 
variations  6v0,  6qj,  and  6qot  and  Eq.  (29)  related  6Gy  and  SGj  with  the  same 
set  of  variations.  These  are  rewritten  as 

26v0  6qt  6q0 

6  Gy  ■  by - +  bj - +•  b0  - -  (50) 

v0  q0-qi  q0-qi 

6q^  ■  -  6q0  +  2m6Gx  —  26Gy  (51) 


By  substituting  Eqs.  (48)  and  (50)  into  (51)  we  have 

26v0  6qi  6q„ 

6qi  -  -  6q0  +  2m (av - +  a* - +  a0 - ] 

v0  qo-qt  q0-qi 

26v0  6qj  6q„ 

-21  by  — —  +  bj  — —  +•  b0  — — ]  (52) 

v0  9o"9i  9o-<li 
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or 


(-2mav+2bv)  26v0  -2a«i+2bi  6qj  -2ma0+2b0  6q0 

0  . . . +  (1  + - ) - +  (1  + - ) - 

qo-qi  vo  qo"qi  q0-<ii  q0-qi  q0"qi 


(53) 


Equations  (49)  and  (53)  are  In  terms  of  the  variations  6v0,  dq*,  and  6q0. 
Sensitivity  coefficients  are  ratio  of  6q0  and  6vQ,  thus  dq*  can  be  eliminated 
by  combining  these  two  equations: 


-2mav+2bv  aj  av  -2maj+2bi  2dv0 

0  -  [( - )<1 - )  4-  (1  +  — )(1  + - )] - 

q0-qi  °x  Gx  q0-qi  vo 


-2ma0+2b0  aj 

+  [(1  + - )(1 - )  + 

qo-qi  g* 


ao  2q0(q0-qi)  -2mai+2bl  6q0 

O  + . rrn - )(1  + - >1 -  (54) 

gx  i+qo  q0-qi  qo-qi 


If  we  define  the  ratio  S*  as 


S* 


«q0/(qo-qi) 

5v0/v0 


Then  from  Eq.  (54)  this  ratio  Is 


(55) 


S*  -  (-2) 


-2mav+2bv  aj  av  -2mai+2b^ 

( - )(1  -  --)  +  (1  +  --)(!  + - ) 

qo-qi  °x  °x  qo-qi 


-2ma0+2b0 


®1  a, 

O  + - )(1  -  --)  +  O  + 

q0-9l  Gx 


o  2q0(q0-q1)  -2ami+2bi 

HI  + - ) 

o,  9o-qi 


(56) 


i+qe2 


The  sensitivity  coefficient  Is  defined  as 

5q0/0+qo2> 

5v0/v0 


(57) 
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which  give  the  relationship 


S 


9o“qi 


(58) 


ITERATION  OF  SOLUTIONS 


The  terrain  slope  m  is  usually  provided  before  the  computation.  Equation 
(27)  shows  that  this  terrain  slope  is  in  terms  of  the  range  drag  function,  Gx 

in  Eq.  (16)  and  the  elevation  drag  function  Gy  in  Eq.  (26).  In  turn,  these 

drag  functions  are  expressed  in  terms  of  the  projectile  starting  slope  qQ  and 
its  Impact  slope  qj.  Computationally  the  solution  of  m  for  given  qQ  and  qj  is 
a  straight  forward  substitution  procedure  using  Eqs.  (27),  (16),  and  (26). 
However,  since  the  terrain  slope  m  is  known  in  advance  physically  and  the 

impact  slope  q^  is  not,  we  compute  instead  the  impact  slope  q*  for  a  given 

starting  slope  qQ  in  Eq.  (27).  Rewriting  this  equation  we  have 


9i 


-  <lo  + 


2 

qo-qi 


qH(q,qo»vo2»c/8>d<i 

"7~-~h 


2m  M  H(q,q0,v02,c/g)dq 

+  2m - /  -  (59) 

q0-qi  q0  1  -  H 


One  can  guess  an  initial  value  for  qj  on  the  right  side  of  the  above  equation 
and  perform  the  integration.  The  value  for  qj  on  the  left  side  can  be  used  as 
the  initial  value  for  q*  on  the  right  side  for  the  second  trial.  This  itera¬ 
tion  procedure  continues  until  the  value  of  qi  converges  to  a  numerical  solu¬ 
tion  as  its  limit. 
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Equation  (59)  may  be  expressed  as 


91 


9o 


2 

+ - 

9o“9i 


•*11 


2a 

+  2a  -  — - —  Ij  i 

9o-9l 


(60) 


where  Ijj  and  Jn  are  given  In  Eqs.  (38)  and  (43).  These  Integrals  will  be 
evaluated  In  the  next  section. 


EVALUATION  OF  INTEGRALS 

If  the  drag  coefficient  cv02/g  Is  small  compared  to  (l+qQ2)  the  denomina¬ 
tor  terms  (1-H)  and  (1-H)2  In  Eqs.  (38)  through  (40)  and  (43)  through  (45)  may 
be  neglected.  One  obtains  from  the  above  equations  the  following 


rqi 

I11  2  *12  -  I  Mq 

9o 

*02  2  <U  “  9o 

Jn  2  Ji2  2  /  1  qMq 

9o 

Jo2  •  -  (qi2-q02) 


Equation  (11)  for  H  may  be  rewritten  as 


H  "  K  TIT”?  ”  Po(9o)l 

where 

K  -  cv02/g 

and  p(q)  Is  given  In  Eq.  (77).  Then  the  Integrals  Ijj  and  Ij2  sre 

*11  “  *12  "  KM 


(61) 

(62) 

(63) 

(64) 

(65) 

(66) 

(67) 
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where 


and 


1 

M  ■  •“"?  [^q«q  "  ^q»q  ^  (qo_<li)Po(9o))  (68) 

l+qQ^  n  H  no 

P  -  /  p(q)dq 

-  -  (l+q2)3/2  +  qin[q  +  (1+q2)1^2]  -  (1+q2)1^2  (69) 

3 


Similarly  the  integrals  Ju  and  Ji2  ate 

Jll  -  Jl2  -  KH  (70) 

where 

N  “  ^Qq*q.  “  Qq*q  +  7  (9o2”<li2)Po(<lo)l  (71) 

l+q0  i  o  2 

and 

Q  ■  J  qp(q)dq 


*  r  q(l+q2)3^2  -  -  q(l+q2)1/2  +  (-  q2  +  “)4n[q  +  (1+q2)1^2]  (72) 

4  8  2  8 


In  the  last  section  the  iteration  equation  (59)  for  the  solution  for  qj 
now  becomes 

qi  *  "  qo  +  2(mGx-Gy)  (73) 


where 

Gx  -  1  -  (qo-qi)’1  KM  (74) 

and 

Gy  -  -  (q0-qi)_l  KN  (75) 

We  are  now  able  to  compute  qj  if  the  quantities  m  and  q0  are  given* 
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The  coefficients  used  in  the  variation  of  range  and  elevation  drag 
functions  in  Eqs.  (32)  through  (37)  give 

av  -  -  (qo^i)"1  KM 

at  -  K{-[p(qi)  -  p0(qo>l  “  (q0"qi)_l  M1 


2qo 


dPo 


ao  -  k{[~2  +  (qo-qi)"l]M  -  ---  (q0-qi)l 

i+q0  dqo 

bv  -  -  (qo-qi)"1  kn 

bi  ■  Ki-qilp(qi)  -  Po(q0>]  ■  (qo~qi)”1  N1 


and 


2q0  I  dp0 

bo  *  k{  [--—j  +  (q0-qi)-1]N  -  -  (q02-qi2)  } 

lo 


l+q02  '  2  --  *  dqQ 


(76) 

(77) 

(78) 

(79) 

(80) 

(81) 


With  qQ  and  q*  known,  we  are  able  to  compute  the  coefficients  of  three  a's  and 
three  b's  as  given  in  the  above  equations.  It  is  now  ready  to  substitute  all 
the  coefficients  into  Eq.  (56)  to  compute  S*.  By  Eq.  (58),  in  turn,  one 
obtains  S,  the  sensitivity  we  are  seeking. 

DETERMINATION  OF  THE  EXTREMAL  VALUE  OF  THE  DAMPING  COEFFICIENT 

The  assumption  that  the  denominator  terms  (1-H)  may  be  dropped  in 
evaluating  the  integrals  depends  on 

H  «  1  (82) 

From  Eq.  (65)  this  is  equivalent  to 


l 

i+q0* 


H(q)  -  k{— j  lp(q)  -  P0(q0)l )  «  1 
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where 

p(q)  ■  {q(l+q2)1^2  +  An(q  +  (1+q2)1/2]}  (84) 

It  Is  noted  that 

H(q0)  -  0  at  q  -  qQ  (85) 

The  extremal  value  of  H(q)  can  be  evaluated  at  q  ■  qj 

|H|  <  |H I extremal  “  lH(q±) I  -  K  ----^  [p(qi)  -  p0(qo)]l  «  1  (86) 

The  above  Is  conservative  In  that  we  use  the  largest  value  of  I H |  for 
determining  the  range  spread  of  K.  For  example  if  qQ  •  1  and  q*  »  -1,  then 

~x  PoUo)  -  O+to  2.414) / 2  -  0.940  (87) 

i+q0 

~--t  p(qt)  -  (-1+An  0.414)/2  -  -0.940  (88) 

i+qQ 


then 


*1-0.940  -  0.940|  «  1 
or 

K  «  0.532 


and  if  qQ  ■  /3,  and  qj  «  -/3,  then 

- - j  p0(q0)  -  1(1 .732) (2)  +  in  3.732]/4  -  1.195 

l+9o 

1 

- -  p(qi)  .  ( (-1 .732) (2)  +  in  0.2681/4  -  -1.195 

i+q0 

then 

*1-1.195  -  1.1951  «  1 


(89) 

(90) 

(91) 

(92) 

(93) 
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K  «  0.418 


(94) 


The  maximum  of  K  we  chose  In  this  paper  is 

max(cv02/g)  ■  max  K  •  0.20  (95) 

which  satisfies  the  above  requirements  for  the  approximation  used  in 
evaluating  the  integrals. 

COMPUTATION  PROCEDURE 

The  computation  procedure  consists  of  two  parts,  the  iteration  procedure 
for  the  Impact  slope  qj  and  the  sensitivity  and  range  computation  after 
obtaining  q^. 

(1)  Iteration  Procedure  -  By  giving  the  launch  slope  qQ  one  can  compute 
such  quantities  as  p0(q0)  from  Eq.  (7),  Pq»q  from  Eq.  (69)  and  Qq«q  from  Eq., 
(72).  The  iteration  procedure  starts  by  assuming  an  initial  Impact  slope  q*. 
Then  we  can  obtain  Pq-q^  and  Qq>q^  from  Eqs.  (69)  and  (72),  respectively.  In 
addition,  we  can  compute  M  and  N  from  Eqs.  (68)  and  (71),  respectively.  By 
assuming  the  value  of  K  (i.e.,  cv02/g)  in  Eq.  (66)  the  integrals  I}}  and 

are  calculated  from  Eqs.  (67)  and  (70),  respectively.  With  these  Integrals 
one  can  find  a  new  value  for  the  Impact  slope  qj  by  using  Eq.  (60).  This 
iteration  procedure  continues  until  the  value  of  qj  converges  to  a  numerical 
solution  as  its  limit. 

(2)  Sensitivity  and  Range  Computation  -  With  the  iterated  solution  of  the 
Impact  slope  qj  known,  one  can  proceed  to  find  the  sensitivity  coefficients 
and  the  range.  We  start  to  find  the  range  drag  function  Gx  and  the  elevation 
drag  function  Gy  from  Eqs.  (74)  and  (75),  respectively.  Next  the  coefficients 
av,  aj,  a0,  bv,  bj ,  and  b0  can  be  evaluated  by  Eqs.  (76)  through  (81). 
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Substituting  these  functions  and  coefficients  into  Eq.  (56)  we  can  get  S*, 
thus  the  sensitivity  coefficient  S  can  be  obtained  readily  from  Eq.  (58). 

The  range  X  can  also  be  computed  by  combining  Eqs.  (13),  (15),  and  (74). 

NUMERICAL  SOLUTIONS 

We  are  using  the  nondimensional  damping  coefficient  K  ■  cv02/g  ■  0.2  for 
illustration  purpose.  The  terrain  slopes  are  assumed  to  be  m  *  0  and  m  - 
0.2679  which  is  corresponding  to  a  terrain  angle  of  0®  and  15°  respectively. 
The  launch  angle  60  are  varied  as  follows: 

for  m  »  0,  7.5°  <  0O  <  82.5°  at  intervals  of  7.5® 
for  tan"1  m  -  15®,  22.5  <  0O  <  75®  at  intervals  of  7.5® 

The  results  are  shown  in  Tables  I  and  II,  and  plotted  in  Figures  2  through  6. 
Figure  2  gives  the  Iterated  Impact  angles  for  various  launch  angles.  The 
effect  of  K  to  the  range  is  fairly  large  as  shown  in  Figures  3  and  4.  How¬ 
ever,  the  effect  of  the  nondimensional  drag  coefficient  K  to  the  sensitivity 
coefficient  is  not  very  large  except  near  the  critical  sensitivity  which  mag¬ 
nitude  is  unbounded  as  shown  in  Figures  5  and  6. 

CONCLUSION 

This  report  is  the  third  of  a  sequence.  The  findings  are  summarized  as 


follows 


From  Che  first  paper 1  we  have: 

(1)  The  error  Increments  in  hitting  a  target  are  derived  and  set  to  zero 
for  a  bull's  eye  landing. 

(2)  The  sensitivity  coefficient  is  defined  as  the  ratio  of  increments  of 
the  initial  elevation  angle  to  the  increments  of  the  natural  logarithm  of  the- 
velocity. 

(3)  Without  air  resistance  or  damping  the  maximum  range  is  found  to  be 
located  at  an  optimal  initial  elevation  angle,  which  is  half  of  the  sum  of  904 
and  the  terrain  angle. 

(4)  Without  damping  the  sensitivity  coefficients  become  infinite  at 
maximum  range  and  zero  at  minimum  range  of  zero. 

(5)  The  sensitivity  coefficient  is  positive  for  high  trajectories  (above 
45°  initial  elevation  angle)  and  negative  for  low  trajectories  (below  45* 
initial  elevation  angle). 

(6)  Due  to  uncertainty  of  muzzle  velocity  deviation  a  gun  may  be  designed 
by  matching  the  interior  ballistics  to  the  exterior  ballistics.  This  is  to 
give  automatically  and  Instantly  a  proportional  deviation  of  its  initial 
elevation  angle  when  the  shell  leaves  the  muzzle  during  firing  for  the  first 
round . 


* Shen ,  C.  N. ,  "On  the  Sensitivity  Coefficient  of  Exterior  Ballistics  and  Its 
Potential  Matching  to  Interior  Ballistics  Sensitivity,"  Proceedings  of  the 
Second  US  Army  Symposium  on  Gun  Dynamics  at  the  Institute  on  Man  and  Science, 
Rensselaervllle,  NY,  September  1978,  sponsored  by  USA  ARRADCOM. 
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The  second  paper^  Investigated  the  problem  when  the  principal  equation  of 
exterior  ballistics  has  a  drag  term  which  is  proportional  to  the  square  of 
velocity.  The  paper  is  to  find  analytically  the  sensitivity  coefficient,  A 
sequence  of  nonlinear  transformation  was  used  before  the  sensitivity  problem 
can  be  studied.  The  variation  of  various  parameters  is  investigated, 
including  the  variation  under  an  integral  sign. 

The  present  report  continues  the  work  of  the  second  and  converts  it  into 
a  different  analytical  form  that  is  suitable  for  numerical  computation.  The 
solution  for  sensitivity  coefficients  is  condensed  into  a  single  formula  with 
the  terms  readily  computable  from  a  set  of  parameters.  An  iteration  procedure 
is  used  in  the  computation.  Some  integrals  are  approximated  in  order  to 
obtain  their  analytical  solution  forms. 

The  effect  of  velocity  square  damping  to  the  range  can  be  large, 
especially  at  the  maximum  range.  However,  the  effect  of  velocity  square 
damping  to  the  sensitivity  coefficient  is  usually  small,  except  near  the 
maximum  range . 


^Shen,  C.  N.,  "Sensitivity  Coefficient  of  Exterior  Ballistics  With  Velocity 
Square  Damping,"  Transactions  of  the  Twenty-Fifth  Conference  of  Army 
Mathematicians,  ARO  Report  80-1,  January  1980,  pp.  267-282. 


21 


TABLE  I.  RANGE  AND  SENSITIVITY  COEFFICIENTS  FOR  FLAT  TERRAIN 


For  m  *  0  cv02 


eo 

9o 

g 

BjCdegree) 

91 

X 

S 

7.5® 

0.13165 

0 

-  7.5000 

-0.131652 

0.258819 

-0.267949 

0.20 

-  7.7808 

-0.136641 

0.249773 

-0.268020 

15.0° 

0.26794 

0 

-15.0000 

-0.267949 

0.500000 

-0.577350 

0.20 

-16.1583 

-0.289738 

0.465538 

-0.578836 

22.5° 

0.41421 

0 

-22.5000 

-0.414213 

0.707106 

-1.000000 

0.20 

-25.0595 

-0.467574 

0.636169 

-1.011472 

30.0° 

0.57735 

0 

-30.0000 

-0.577350 

0.866025 

-1.732050 

0.20 

-34.2094 

-0.679841 

0.755507 

-1.804402 

37.5° 

0.76732 

0 

-37.5000 

-0.767326 

0.965925 

-3.732050 

0.20 

-43.1708 

-0.938106 

0,821550 

-4.422821 

45.0° 

1.00000 

0 

-45.0000 

-1.000000 

1.000000 

00 

0.20 

-51.4995 

-1.257153 

0.835103 

14.3160 

52.5® 

1.30322 

0 

-52.5000 

-1.303225 

0.965925 

3.732050 

0.20 

-58.9394 

-1.660305 

0.798363 

2.645294 

60.0® 

1.73205 

0 

-60.0000 

-1.732105 

0.866025 

1.732050 

0.20 

-65.5374 

-2.198108 

0.714006 

1.342668 

67.5® 

1.73205 

0 

-67.5000 

-2.414213 

0.707106 

1.000000 

0.10 

-69.0956 

-2.618147 

0.647751 

0.915074 

0.20 

-71.6039 

-3.006806 

0.585294 

0.796916 

75.0® 

3.73205 

0 

-75.0000 

-3.732050 

0.500000 

0.577350 

0.10 

-76.0085 

-4.013342 

0.459463 

0.529802 

0.20 

-77.5372 

-4.524636 

0.417197 

0.465626 

82.5® 

7.59575 

0 

-82.5000 

-7.595754 

0.258819 

0.267949 

0.10 

-82.9614 

-8.099269- 

0.238697 

0.245899 

0.20 

-83.6382 

-8.969194 

0.217891 

0.217213 
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TABLE  II.  RANGE 


SENSITIVITY  COEFFICIENTS  FOR  15*  TERRAIN 


AND 


For  a 
60 

-  +0.2679 
9o 

cvo2_ 

8 

®l(degr««) 

<U 

X 

S 

22.5* 

0.41421 

0 

6.932362 

0.121586 

0.249772 

-0.278596 

0.10 

6.783131 

0.118944 

0.245353 

-0.278737 

0.20 

6.619309 

0.116045 

0.240893 

-0.278891 

30.0* 

0.57735 

0 

-  2.379286 

-0.041550 

0.464175 

-0.634112 

0.10 

-  3.002055 

-0.052443 

0.448368 

-0.636174 

0.20 

-  3.751089 

-0.065562 

0.432262 

-0.638617 

37.5* 

0.76732 

0 

-13.03583 

-0.231526 

0.628688 

-1.214643 

0.10 

-14.38786 

-0.256530 

0.598311 

-1.230168 

0.20 

-16.14129 

-0.289416 

0.567092 

-1.250663 

1.00000 

0 

-24.90070 

-0.464200 

0.732100 

-2.732736 

0.10 

-26.98236 

-0.509137 

0.688258 

-2.877936 

0.20 

-29.80982 

-0.572933 

0.642910 

-3.107193 

52.5* 

1.30322 

0 

-37.50354 

-0.767425 

0.767364 

«0 

0.10 

-39.96784 

-0.838143 

0.714885 

32.535621 

0.20 

-43.37056 

-0.944680 

0.660410 

13.127275 

60.0* 

1.73205 

0 

-50.10622 

-1.196250 

0.732075 

2.731708 

0.10 

-52.42571 

-1.299732 

0.678076 

2.405832 

0.20 

-55.62687 

-1.461933 

0.622002 

2.031230 

67.5* 

2.41421 

0 

-61.97076 

-1.878413 

0.628640 

1.214358 

0.10 

-63.73921 

-2.026837 

0.580688 

1.105513 

0.20 

-66.17359 

-2.264476 

0.531007 

0.968816 

75.0* 

3.73205 

0 

-72.62684 

-3.196250 

0.464108 

0.633962 

0.10 

-73.71086 

-3.422142 

0.428564 

0.581436 

0.20 

-75.20812 

-3.787022 

0.391874 

0.514330 

82.5* 

7.59575 

0 

-81.93802 

-7.059954 

0.249690 

0.278489 

0.10 

-82.41359 

-7.508239, 

0.230773 

0.255652 

0.20 

-83.41359 

-8.238227 

0.211307 

0.226622 
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APPENDIX  A 


DYNAMICAL  EQUATIONS  FOR  TRAJECTORIES  AND  VARIABLE  TRANSFORMATIONS 

For  a  constant  mass  travelling  in  a  vertical  plane  with  no  life  and 
applied  thrust,  but  having  drag  and  velocity  vectors  contained  in  the  plane  of 
symmetry  as  shown  in  Figure  1,  the  dynamical  equations  of  motion  are:1 

dx  ,  v 

—  -  v  cos  8  ■  0  (Al)  , 

dy 

—  -  v  sin  8  ■  0  (A2) 

dt 

d6 

m(g  cos  8  +  v  —  )  ■  0  (A3) 

dt 


d2x  D  cos  0 


It  is  noticed  that  deviations  due  to  anomalies  in  the  azimuth  direction  are 
not  considered  here. 

By  differentiating  Eq.  (Al)  with  respect  to  t  one  obtains 


dzx  d 

— r  ■  —  (v  cos  8) 
dt2  dt 


Substituting  into  Eq.  (A4)  we  have 


d  D  cos  6 

—  (v  cos  0)  ■  - - — 

dt  m 
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Solving  for  dd/dt  In  Eq.  (A3)  one  obtains 


d8  -g  cos  8 
dt  v 


(A7) 


Equation  (A7)  indicates  that  the  differential  equations  can  be  transformed 
from  the  time  domain  in  t  to  the  angle  domain  in  8.  Equations  (A6),  (Al),  and 
(A2)  are  divided  by  Eq.  (A7)  in  achieving  this  transformation  as 


d(v  cos  8)  Dv 
dd  mg 


(A8) 


dx  v2 

de  “  “  g’ 


dy  v2 

—  ■  .  —  tan  8 
d8  g 


(A9) 


(A10) 


Equation  (A8)  is  called  the  prlnicpal  equation  of  exterior  ballistics.2  It 
can  be  integrated  if  the  drag  D  is  a  known  function  of  velocity  v. 

The  head  wind  drag  is  a  velocity  square  damping  term  given  as 

D  •  mcv2  (All) 

where 

c  -  cw(^  d2)(p/2)  (A12) 

4 

cw  -  the  dimensionless  resistant  coefficient 
d  -  the  diameter  of  projectile 
and  p  ■  the  air  density. 
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Thus  Che  principal  equation  of  exterior  ballistics  (Eq.  (A8))  becomes 


d  cv*  , 

—  (v  cos  0)  ■ -  (A13) 

d6  g 

A  further  transformation  of  the  dependent  variable  in  necessary  by  letting 

u  ■  v  cos  0  (A14)  ■ 

where  u  is  the  horizontal  component  of  the  projectile  velocity.  Then  the 
dynamical  Eqs.  (A13),  (A9),  and  (A10)  become 

du  c  .  . 

—  ■  -  u3  sec3  0  (A15) 

dfl  g 

dx  u2  „ 

—  -  -  —  sec2  0  (A16) 

dO  g 

dy  u2  , 

—  »  -  —  sec2  0  tan  0  (A17) 

d0  g 

To  simplify  further  the  form  of  the  dynamical  equations  another  transformation 
of  the  independent  variable  is  made  by  letting 

--  -  sec2  0  -  i+q2  (A18) 

«  A  * 


\ 


A-3 


fig.  1.  Forces,  Slopes,  end  Initial  end  Final 
Par  maters  for  a  Trajectory 
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APPENDIX  B 


VARIATION  OF  THE  DRAG  FUNCTIONS 

As  given  by  Eq.  (54)  in  Che  previous  paper?  the  variation  of  the  range 


drag  function  dG*  is 


5G* 


Jl2  ?fiv0 
<Jo-qi  vo 


+ 


Vqi  in  ^  «qi 
i-Hq-qi  q0-qi  q0-qt 


+  [ 


2qoIl2  Ill  c  vo2 

i+qQ2  qo-qi  a  i+q<>2 


dPoI02  ^  «q0 
dq0  q0-qi 


(Bl) 


It  is  noted  that  the  difference  of  the  end  slope  is  not  zero,  l.e.,  qQ-qi  *  0. 
Therefore,  the  problem  does  not  become  singular.  He  have  expressed  the 
variation  6G*  in  terms  of  the  variational  parameters. 

As  given  by  Eq.  (77)  in  the  previous  paper?  the  variation  of  the 
elevation  drag  function  dGy  is 


Jl2  2«v0 

«Gy  "  [ - ] - 

qo-qt  *o 


qiHq-qi  Jn  ^»i 
l_Hq-qi  q0-qi  q0-qi 


♦  [ 


?qod12  J11  c  v02  dpoJ02 

i+q02  qo-qi  s  i+q02  dPo 


dqo 

qo-qi 
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figure  5.  Sensitivity  Coefficient  vs.  Launch  Angle  for 

Horizontal  Target  and  Different  Damping 

Coefficients. 
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Figure  6. 


Sensitivity  Coefficient  vs.  Launch  Angle  for 
Upgrade  Targets  and  Different  Damping  Coefficients. 
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